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SUMMARY 


The  abridged  Nicholson  asympototic  formulas  for 
the  modified  Bessel  functions  of  the  first  and  second 
kinds  are  herein  replaced  by  more  accurate  asympototic 
formulas.  These  are  derived  and  then  applied  to  the 
evaluation  of  illustrative  numerical  examples  of  slowly 
convergent  Kapteyn  series  involved  in  the  calculation 
of  propeller  induction  factors.  The  resulting  values 
of  the  induction  factors  are  then  compared  with  those 
obtained  by  means  of  the  Nicholson  formulas. 


THE  CALCULATION  OF  PROPELLER  INDUCTION  FACTORS 


Introduction 

In  determining  the  axial  and  tangential  components  of  the 

velocity  induced  at  a  lifting  line  by  the  vortex  system  of  a 

(1)* 

propeller,  Lerbs  found  it  convenient  to  use  the  induction 

(2) 

factors  introduced  by  Kawada. 

The  computation  of  these  factors  involves  principally 

the  numerical  evaluation  of  two  infinite  series  of  products 

of  modified  Bessel  functions  and  their  derivatives,  each 

depending  on  three  parameters.  These  series  are  of  the 

Kapteyn  type,  involving  summation  over  both  the  arguments 

and  the  orders  of  the  constituent  Bessel  functions,  and 

consequently  offer  formidable  computational  difficulties 

for  certain  combinations  of  the  parameters  appearing  therein. 

The  numerical  results  obtained  by  Lerbs  were  derived 

through  the  use  of  the  leading  term  of  each  of  Nicholson’s 

(3) 

asymptotic  formulas  for  the  modified  Bessel  functions 

of  the  first  and  second  kind,  namely,  and  . 

~  (*+) 

These  formulas  have  been  aptly  described  by  Lehmer  as 

expressed  "in  terms  of  an  infinite  differential  operator 

with  undetermined  coefficients  operating  on  a  cumbersome 

functiono"  Accordingly,  Lehmer  (loc.  cit.)  derived  an 

♦Numbers  refer  to  references  appearing  bn  page  17  of  this 
report.’ 


alternative  asymptotic  formula  for  which  is  con¬ 

venient  to  apply  to  calculations  requiring  higher  accuracy 
than  that  afforded  by  the  corresponding  abridged  Nicholson 
formula.  The  -formula  obtained  by  Lehmer  suggested  to  the 
present  writer  an  analogous  formula  for  Kjny.), 

which  has 

proved  equally  satisfactory  in  numerical  applications. 

These  improved  formulas  for  XJtv#)  and  ty*#)  permit 
a  more  accurate  approximation  of  the  sums  of  the  Kapteyn 
series  under  consideration. 

It  is  the  purpose  of  this  note  to  derive  the  more 
accurate  formulas, to  illustrate  their  use  by  appropriate 
examples,  and  to  give  a  close  estimate  of  the  error  inherent 
in  the  Nicholson  formulas  as  used  by  Lerbs. 

As  a  result  of  this  study  the  conclusion  is  reached  that 
Lerbs's  numerical  values  for  the  induction  factors  are 


generally  correct  to  within  two  units  in  the  third  decimal 


place,  but  occasionally  may  be  in  error  by  five  units  in 
that  place.  Consequently,  recalculation  of  his  data  seems 
to  be  unnecessary  for  all  practical  applications  considered 


at  this  time. 

Derivation  of  Formulas 

Meissel’s  first  extension  of  Carlini’s  formula  for 


form 


(5) 


as  reproduced  by  Watson 


can  be  written  in  the 
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77  „  \  %  e  usfyfvi—V-n) 


(1) 


where 


V  =  *»;»* 


— f/. _ ikzimizj^iizmi^  (s) 

I  in.’- (I  ~f-)3  irjioJi  (\ -p*}1"* 


We  set  ^  =s  and  use  the  relation 

then  equations  (1)  and  (2)  yield  the  asymptotic  formula: 


yu 


where 

-VftvJ 


TV 


-3>i —  4  ^  i  5’75x-3frh?+i'SfZ±'+  /£ 


+ 


/3*-  —  ~f~3.S8%? — j 


i  as -n.*  (i +  *.*-)  ** 


(4) 


(3) 


C4) 


(5) 


An  examination  of  Lehmer 's  derivation  reveals  that 

formula  (4)  is  valid  for  all  real  values  of  jC  . 

The  analogous  formula  for  K is  most  expeditiously 

^  .  ,ro 

derived  from  the  expansion  obtained  by  Meissel  for  H  (y^aju^^) 
in  conjunction  with  the  fundamental  relation 


3 


(6) 


which  is  given  by  Watson  (loc.  cit.,  p.  78) 
Meissel 's  formula  is 


e. 

* 


(7) 


where 

p  =  +^/s; 

-  (SA  +*S*./>*£p 


<Q  =  7i,(tL~&-p)  -  (&  +  3#J-6) 

'  '  3L+V*  ' 

— 3  (/£  —iSia.A*&p  -‘3(»slAufjp  —3ySj»t <£j$) 

S7«0TV  (9j 


We  set  sec  ^  =  i  #  in  equations  (7)  and  X8),  and 
then  use  equation  (6)  to  obtain  the  desired  asymptotic 


formula,  namely 


,  do) 


(11) 


where  we  have  followed  Lehmer's  notation  in  setting 


=  &(***-*) 
Vj  (*■*)  —  7 £-(*?-**} 


-i  (13*?-  J23a.*f  -f-J-tt  *?-3a*a,\ 
US  V  ' 


(4) 

Lehmer  has  published  four  additional  polynomials  of  this 


family,  and  has  given  formulas  for  the  determination  of 
additional  ones,  if  they  are  desired. 

The  similarity  of  formulas  Of)  and  (10)  is  more  clearly 


brought  out  if  it  is  noted  that 


in  formula  Of). 


(12) 


The  asymptotic  formulas  just  derived  can  be  used  to 


yield  similar  formulas  for  the  logarithmic  derivatives  of 
the  modified  Bessel  functions. 


We  infef  the  existence  of  asymptotic  expansions  for 
I  (rt/)C)  and  from  the  fact  that  these  deriva- 

W 


tives  are  expressible  as  linear  combinations  of  Bessel 
functions,  as  exhibited  by  the  formulas: 


(13) 
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Kjrt)  =  -klKJV  + 

(6) 

According  to  a  theorem  given  by  Knopp  we  can  deduce  the 
series  for  the  logarithmic  derivatives  of  lb*#.)  and 
K  (‘A^J  by  means  of  term-by-term  differentiation  of  the 

Vw 

logarithms  of  the  series  in  equations  (*0  and  (10),  respec¬ 
tively. 


(1*0 


The  desired  results  are 

~  k±*£-  t.  + 

ru  *  An^l4*.7 


iWl 

ru  1  ru 


f)' 


~  2r}(lZ^T~  +  tki-rSO  +  tf^ — 

Km/kw  ~  -a±i‘-  •  -  £te&r 


*.(4  -  IQ'L'+tf)  '  *-( W—Stoff+y#  it-Asd) _ (16) 


t^O-h^r 


We  are  now  in  a  position  to  derive  asymptotic  formulas 
for  the  terms  of  the  two  Kapteyn  series  under  consideration. 
Explicitly,  these  series  are 
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(17) 


f.  *  £  ^ 


In  these  stuns  g  represents  a  positive  integer,  usually  3,  k 
5,  or  6  (corresponding  to  the  number  of  blades  in  a  given 
propeller),-^  is  a  positive  number,  and  ^  is  a  second 
positive  number  such  that  ^  in  equation  (17)  and 

>*  -^0  in  equation  (18). 

Equations  (lV) ,  (10),  and  (16)  imply  the  asymptotic 
relation. 

-  ]>T 1> ~  -  ^  (-3$  {’*&-  jgtri 


A  C 
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l  + 


'f  vfy 

+  4^+  +" 
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(19) 


Similarly,  equations  (4),  (10),  and  (15)  imply  the 
asymptotic  relation 
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^  ^T- -r 


jwf 


T~ 


m*l 


S'^S'^’t)  -  -h  irO  ^  7 

wtyyr  fyi+ffiri 


(20) 


The  Nicholson  approximations,  as  used  by  Lerbs,  lead  to 


the  approximate  results 

/  y  -y^/  £ 

nsl(nfif)  KWfcJ/v  _  — ! — _ /ii&_)  H  j  I  +  — & — ™-  1  (21) 
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The  members  of  formulas  (21)  and  (22)  can  be  summed  in 
closed  form  to  yield  the  following  formulas  used  by  Lerbs 


in  his  calculations 


The  inherent  errors  in  these  formulas  are  closely  approxi¬ 
mated  by  the  following  additive  corrections: 


_  _J _ 

for  f; 


I  + 


-J-) 

U-// 


(27) 


for  F,  . 

X 

'J’hese  corrections  can  be  deduced  by  appending  to  the  Nicholson 
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formulas,  (21)  and  (22),  the  first  two  terms  of  the  Maclaurin 
development  of  the  final  exponential  terms  appearing  in 
equations  (19)  and  (20).  *411  terms  involving 

to  higher  powers  than  the  first  are  discarded  prior  to 
summation.  Inclusion  of  such  terms  leads  to  the  summation  of 

n  /  tt, 

series  whose  general  terms  are  of  the  form  U,  fyy,  y  where 
r  £  2,  and  such  results  are  expressible  in  terms  of  Spence's 
integrals,  which  are  not  tabulated  except  for  r  —  2. 

Nicholson's  formulas  when  modified  by  the  respective 
corrections  shown  in  (26)  and  (27)  can  be  written  in  the 
forms : 


F 

i 


(0Ai^ 


_i_r  q 

(i + 


(28) 


(29) 
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Illustrative  Calculations 

The  evaluation  of  and  when  ^  and  ^  are  nearly 
equal  is  complicated  by  the  corresponding  slow  convergence 
of  the  respective  Kapteyn  series.  Indeed,  when  ^  and  ^ 
are  equal  these  series  diverge. 


We  illustrate  this  difficulty  and  the  manner  in  which 
it  is  overcome  by  computing  rj  corresponding  to  the  following 
data:  g  s  3,  y  a  10/  11,  ^  as  1. 

In  this  example  we  can  calculate  and 

(7) 

for  n  *  1,  2,  ...,  6  by  means  of  the  extensive  tables 
published  for  the  British  Association  for  the  Advancement 

of  Science.  The  derivative  K.  (^y\J  was  calculated  by 

syu 

means  of  eq  (14).  Full  accuracy  in  interpolation  was 
attained  by  using  Everett's  interpolation  formula  with 
second  central  differences. 

We  tabulate  next  these  values  of  yu  X  j 


3>v 


5  >v 


and  for  comparison  show  beside  them  the  corresponding  data 


computed  by  means  of  the  formula 


11 


which  constitutes  the  approximation  underlying  formula  (28). 

It  is  informative  to  include  in  the  table  also  the  appropriate 
values  computed  with  the  aid  of  formula  (21). 


Table  of  values  of 


n 

Accurate  value 

Approx,  derived 
from  (30) 

Approx,  derived 
from  (21) 

1 

-0.121165 

-0.121349 

-0.121603 

2 

-0.079521 

-0.079543 

-0.079628 

3 

-0,053079 

-0.053084 

-0.053123 

4 

-0.0355^6 

-O.93559O 

-0.035609 

5 

-0.023904 

-0.023905 

-0.023915 

6 

-0.0l607i 

-0. 016071 

-0.016077 

These  data  are  sufficient  to  permit  an  accurate  evalu¬ 
ation  of  Fj  to  5  decimal  places.  This  is  accomplished 
by  the  following  procedure.  Formula  (28)  is  employed  to 
sum  the  approximating  series  whose  initial  terms  occupy 
the  second  column  in  the  preceding  table.  The  numerical 


result  obtained  is  F  s’  -0.362610.  The  difference  between 


corresponding  entries  in  the  first  two  columns  of  the  table 


form  a  rapidly  converging  series,  namely  0. 000184+  0.000022 


+  0.000005  +  0.000003  +  0.000001+  ...,  whose  sum,  0.000215, 
when  added  to  the  foregoing  estimate  of  yields  the 
accurate  value  F,  =  -0.36239* 

We  compare  this  result  with  that  obtained  by  formula  (24  )^ 


12 


which  yields  the  approximation  Ft  —  -0.36303*  The  latter 
is  seen  to  be  numerically  too  large  by  about  0.18  per  cent. 
The  induction  factors  for  the  internal  field  are  then 


evaluated  by  the  formulas: 


(3D 


(32  ;> 


Substitution  of  the  correct  value  of  F  gives  the  results 

X  .  0. 8657  and  X  s:  0.6523.  whereas  the  value  derived 

eu 

by  formula  (24)  yields  A.  .  —  0.8668  and  a.  at  0.6535, 

jtc 

which  are  therefore  in  error  by  only  0.13  and  0.18  per  cent, 


respectively. 

We  next  illustrate  the  evaluation  of  F  when  the 

•3L 

associated  Kapteyn  series  is  slowly  convergent.  The  specific 
values  of  the  parameters  are:  g  =  3,  y  =  «^/0.9 75,  and 
^  ==•  to*  *  .  The  calculation  of  ^  then  involves 

and  I(WFJ  corresponding 


the  evaluation  of 


K 

3~n 


to  successive  small  integer  values  of  n.  As  in  the  preceding 
calculation,  we  can  use  the  British  Association  tables  to 


compute  these  functional  values  for  n  ranging  from  .1  to  6, 
and  concurrently  evaluate  these  numbers  by  formula  (22)  and  by 
the  following,  more  accurate  formula  which  underlies 


13 


formula  (29). 


By  these  methods  the  following  tabular  data  were  found. 

Table  of  Values  of  K  ,)  J 

3rb  0.97,5V  3-^ 


n  Accurate  value  Approx,  derived  Approx,  derived 


from  (33) 

from  (22) 

1 

0.253838 

0.253835 

0.254001 

2 

0.236885 

0.236886 

0.236962 

3 

0.2182  86 

0.218286 

0.218333 

4 

0.200520 

0.200520 

0.200552 

5 

0.183971 

0.183971 

0.183995 

6 

0.168684 

0.168684 

0.16 87 02 

Evaluation  of  formula 

■  (29)  yields  the 

approximation 

^  3.102481, 

whereas 

the  correct  value 

appears  to  be 

3.: 

102483,  after  making  a 

comparison  of  corresponding  entries 

in 

the  first  two 

columns 

of  the  last  table 

.  These  results  are 

to  be  compared  x^ith  the  approximation  3.102931  obtained 

by  use  of  formula  (25),  which  consequently  seen  to  be  too 
large  by  only  0.014  percent. 

14 


The  associated  induction  factors  for  the  external  field 


are  then  calcilated  by  means  of  the  formulas! 


i  s  A 

OLl. 


f)  F* 


*tju° 


(3^) 


(35) 


The  accurate  value  of  jF  yields  the  results  0.4773 

and  X.  s s  0.8810,  whereas  the  value  of  F  obtained  by  use 

JtjL 

of  formula  (25)  yields  the  approximations  i  —  0.4774  and 

ojl 

l  «  0.8812. 

It  should  be  observed  that  formulas  (28)  and  (29)  apply 
with  equal  facility  to  the  evaluation  of  I6,  and  when 
^  and  ^  are  not  nearly  ecual.  For  example,  when  3, 

s  i^,  and  i.  formula  (28)  yields  the  approximation 
^  m  -0.016063  instead  of  the  accurate  value  Fj  =  -0.016035, 
obtained  with  the  aid  of  the  British  Association  tables. 


On  the  other  hand,  the  abridged  Nicholson  formula  [ equation 
(24 )]  gives  the  approximation  -0.016321,  which  corresponds 

to  1*6469  and  0*2938,  as  contrasted  with  the 

accurate  values  1.6443  and  0.2886,  respectively.  It  is 

interesting  to  observe  that  the  relative  error  an  i  .  if 

cu 

only  0.16  percent,  whereas  the  relative  error  in  ^  is  1.8 
percent.  The  respective  absolute  errors  are  seen  to  be 
2.6  x  10  ^  and  5*2  X  10 

A  similar  calculation  of  corresponding  to  the  para- 

15 


metric  values  g*  3>  y  *  2,  and  1  gives  the  approxi¬ 
mation  0,0032768  when  the  Nicholson  formula  [equation 

(25)J  is  used.  The  accurate  value,  0.0032636,  is  the  sum  of 
the  rapidly  convergent  series:  0.0031?  752  -f  0.00008  391 
+  0.00000217  ~h  ...,  derived  from  data  supplied  by  the 
British  Association  tables.  Formula  (29)  gives  the  result 
~  0.0032510,  which  when  substituted  in  formulas  (3*0 

and  (35)  gives  i  -  0.0585  and  X  m  1.5293. 

ax  Jtx 

For  purposes  of  comparison  we  note  that  the  accurate  values 

are  x  s  0.0587  and  i  m  1.529**»  while  the  values  derived 

ax  tx 

from  the  Nicholson  approximation  to  F  are  respectively 

J4. 

0.0590  and  I.5295,  thus  corresponding  to  errors  of  3-xlO 
-4 

and  1x10  ,  respectively. 

Conclusions 

The  preceding  calculations  indicate  the  reliability  to 
within  5  x  10“3  of  approximations  obtained  for  propeller 
induction  factors  by  means  of  the  Nicholson  formulas  as 
used  by  Lerbs.  The  inherent  errors  can  be  reduced  by  an 
order  of  magnitude  if  formulas  (28)  and  (29)  are  used 
instead.  Finally,  these  small  residual  errors  can  be 
removed  (to  at  least  five  decimal  places)  if  the  first  one 
or  two  terms  of  the  series  defining  F|  and  are  calculated 
independently  by  tables  or  by  suitable  power  series  and  then 
the  corresponding  corrections  are  made  in  the  results  obtained 
by  the  approximating  formulas. 
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